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SUMMARY 


An  analytic  method  for  calculating  the  standard  deviation  of  cloud 
material  vertical  distribution  (sigma  z)  for  a  surface  source  over  short  to 
moderate  travel  distances  has  been  developed.  The  surface  source  sigma  z  is 
described  as  the  product  of  travel  distance  and  the  vertical  component  of 
atmospheric  turbulence,  modified  by  a  correction  factor  (i).  The  resultant 
sigma  z  equation  permits  calculation  of  cloud  vertical  growth,  from  meteor¬ 
ological  profile  data,  over  a  continuum  of  stability  and  travel  distance. 

The  correction  factor,  which  is  a  measure  of  departure  from  rectilinear 
cloud  growth,  assumes  the  form  of  an  arctangent  equation.  The  general  form 
of  the  arctangent  equation  is  1  =  (-ATAN(a*Bl!  +  b)  +  c)/d.  Stability  is 
represented  by  bulk  Richardson  number  (BU),  and  coefficient  a  is  a  constant. 
Variables  b,  c,  and  d  are  functions  of  travel  distance.  Coefficients  a,  b,  c, 
and  d  were  chosen  to  provide  a  best  fit  field  of  i  values  which  give  sigmas 
satisfying  Gaussian  equations  for  Project  Prairie  Grass  time  averaged  cross- 
wind  integrated  concentration  data. 

The  vertical  component  of  atmospheric  turbulence  for  lapse  conditions  is 
shown  to  be  a  function  of  a  stability  parameter,  z/L,  which  is  obtained  from 
meteorological  profile  data.  During  inversion  conditions,  turbulence  is 
described  by  means  of  a  regression  equation  using  bulk  Richardson  number. 

The  sigma  z  computation  procedure  has  been  successfully  used  on 
independent  sets  of  diffusion  data  collected  at  the  U.  S.  Army  Dugway 
Proving  Ground,  Utah,  and  the  National  Reactor  Testing  Station  (NRTS) ,  Idaho. 
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I.  INTRODUCTION 


Gaussian  cloud  models  used  for  the  prediction  of  dosages  or  concentra¬ 
tions  downwind  from  a  surface  source  are  critically  dependent  on  the 
determination  of  a  representative  transport  wind  speed  (u)  and  on  an  adequate 
estimation  of  cloud  dimensions.  Cloud  dimensions  in  Gaussian  models  are 
represented  by  cloud  material  standard  deviations  (sigmas).  Gaussian  models 
are  useful  for  prediction  If  the  u  and  sigmas  can  be  estimated  from  meteor¬ 
ological  profile  data.  Transport  wind  speed  can  be  readily  estimated  from 
measured  wind  profile  data  If  the  cloud  vertical  dimensions  are  known,  but 
present  procedures  for  estimating  cloud  dimensions  from  meteorological 
profile  data  are  not  adequate.  For  a  surface  source,  calculation  of  cloud 
material  vertical  standard  deviation  (sigma  z)  has  been  particularly 
difficult  due  to  the  steep  vertical  gradients  of  temperature,  wind  speed, 
and  turbulence  which  frequently  occur  near  the  Earth's  surface. 

This  report  descrit  -  the  development  of  an  analytic  method  for  using 
estimates  obtained  from  meteorological  profile  data  of  the  vertical  component 
of  atmospheric  turbulence  to  calculate  sigma  z.  Sigma  z  is  the  product  of 
travel  distance  and  the  vertical  turbulence  estimate,  modified  by  a  correction 
factor  (1).  This  correction  factor  accomodates  departures  from  Gaussian 
assumptions  and  turbulence  measurement  deficiencies  over  a  continuum  of 
stability  and  downwind  travel  distance. 

II.  DEFINITIONS 


The  following  symbols  for  standard  variables  and  dimensionless 
quantities  are  defined  for  use  In  the  discussion  to  follow: 


g  =  9.8  m/sec2 
k  =  0.41 
z,  meters 
z0,  meters 
u,  w,  m/sec 
p,  g/m3 

T,  degrees  Kelvin 
9,  degrees  Kelvin 
u*,  m/sec 
0*  =  w?/ku* 
z/L  =  kgwVz/eui 


acceleration  due  to  gravity 
von  Karman's  constant  (dimensionless) 
height  of  measurement  above  the  surface 
roughness  length 

longitudinal  and  vertical  wind  components 
air  density 

absolute  air  temperature 
potential  temperature 
friction  velocity 

scaling  temperature  (degrees  Kelvin) 
dimensionless  height,  a  stability  parameter 
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Ri  =  (g/0)(39/9z)/(3u/32)2 


Richardson  number,  dimensionless 
BU  =  100z2(g/e)(39/3z)/u2  bulk  Richardson  number,  dimensionless 

<|>m  =  (kz/u*)(3u/3z)  dimensionless  wind  shear 

<f>h  =  ( z/e^ ) ( 36/3z )  dimensionless  temperature  gradient 

An  overbar  Indicates  a  time  averaged  quantity  and  a  prime  denotes 
fluctuations  about  the  average.  The  Ri  and  BU  are  stability  parameters 
which  apply  at  the  geometric  mean  height  between  top  and  bottom  of  the 
layer  considered. 

III.  THE  GAUSSIAN  EQUATION,  TRANSPORT  WIND,  AND  SIGMA  Z 

The  analytic  method  was  developed  by  satisfying  the  Gaussian  crosswind 
integrated  dosage  equation  transport  wind  and  sigma  z  requirements  using 
meteorological  profile  data.  Equation  (1)  is  a  Gaussian  equation  for  the 
crosswind  integrated  dosage  of  non-depositing  tracer  material  measured  at 
height  z.  Total  reflection  from  the  surface  is  assumed  using  the  method  of 
images  presented  by  Sutton  (1953). 


Q[EXP(-!s((z-h)/o.)2)+EXP(-J5((z+h)/a  J2)] 

CWID  =  - : - - - - -  (1) 

(2.)%zu 


where: 

CWID  =  crosswind  integrated  dosage  (mg*sec*m-3) 

Q  *  source  strength  (milligrams) 

h  =  release  height  (meters) 

o ^  =  sigma  z 

Strong  gradients  of  wind  speed  frequently  exist  near  the  around. 
However,  if  the  distribution  of  the  wind  speed-dosage  product  (D^u-jdz)  is 
normal  with  height,  then  equation  (2)  defines  a  weighted  mean  wind  speed 
which  can  be  used  as  a  transport  wind  without  violating  the  assumptions 
needed  to  Integrate  the  Gaussian  equations.  Estimates  of  transport  wind 


u  =  /Dudz//Ddz 


(2) 


speed  from  meteorological  profile  data  were  obtained  using  equation  (3) 
developed  by  H.  E.  Cramer  e_t  (1972).  Equation  (3)  uses  the  2-meter 
wind  speed  for  ur  during  daytime  (lapse)  conditions,  and  the  1-meter  wind 
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speed  for  nighttime  (surface  based  Inversion)  conditions.  Variable  zr  is  the 
height  of  the  ur  measurement  and  2.15oz  is  the  "top"  of  the  cloud,  orthe 
height  at  which  cloud  dosage  is  reduced  to  one  tenth  of  its  centroid  value. 
Exponent  p  is  the  wind  speed  profile  exponent  taken  from  zr  to  16  meters. 


uJ(2.15o7)(l+p)  -  z  ^1+P)l 


u  = 


(2.15az-zr)zP(l+p) 


(3) 


With  u  specified  by  equation  (3)  and  the  other  variables  measured, 
equation  (1)  can  be  used  to  solve  for  sigma  z  by  means  of  a  Newton-Raphson 
Iterative  method.  But  before  this  is  done,  a  careful  examination  of  contin¬ 
uity  considerations  is  required.  Because  diffusion  problems  involve  the 
release  of  a  gas  or  a  large  number  of  minute  particles,  the  source  strength 
is  often  imprecisely  known.  A  further  continuity  consideration  is  the 
deposition  rate  or  physical  decay  for  materials  which  have  limited  residence 
time  in  the  atmosphere.  With  a  passive  gas  or  minute  particles  traveling 
over  short  distances,  it  can  be  assumed  that  the  total  amount  of  material 
released  into  the  atmosphere  remains  airborne.  For  other  materials, 
diffusion  equations  must  be  modified  to  include  decay  or  deposition  terms. 
Sampler  collection  efficiencies  and  losses  due  to  data  handling  and  assay 
techniques  are  also  continuity  considerations.  Unfortunately,  these  factors 
are  often  not  carefully  specified  for  diffusion  experiments.  Because  of 
the  difficulties  in  satisfying  continuity  considerations,  the  amount  of 
hi  oh  quality  data  providing  information  required  for  the  solution  of  equation 
(l)  is  limited.  Project  Prairie  Grass  data  (Haugen,  1959)  collected  in  1957 
at  O'Neill,  Nebraska,  were  chosen  for  equation  (1)  sigma  z  solutions  because 
of  the  quality  and  quantity  of  the  data. 

IV.  PROJECT  PRAIRIE  GRASS 

Project  Prairie  Grass  was  a  field  diffusion  program  designed  to 
investigate  diffusion  from  a  continuous  point  source  over  flat  terrain. 

For  a  5  to  6  centimeter  grass  cover,  Pasquill  (1978)  estimated  a  surface 
roughness  of  0.8  centimeters.  The  source  consisted  of  a  nozzle  emitting 
sulphur  dioxide  (S02)  at  a  height  of  either  0.5  or  1.5  meters  above  ground 
level.  Source  duration  was  ten  minutes  for  each  diffusion  trial.  Samplers 
set  out  in  arcs  50,  100,  200,  400,  and  800  meters  from  the  source  provided 
S02  concentration  data  measured  at  a  height  of  1.5  meters.  The  100-meter 
arc  also  included  a  vertical  sampler  array  for  gathering  cloud  vertical 
concentration  measurements. 

The  question  of  continuity  for  Project  Prairie  Grass  data  was  treated 
through  careful  source  strength  measurements  and  by  assuming  that  the 
diffusing  S02  remained  airborne  as  It  moved  across  the  sampling  grid. 

Barad  and  Haugen  (1959)  and  Elliott  (1961)  considered  the  question  of 
deposition  on  vegetation  and  concluded  that  virtually  all  of  the  released  S02 
remained  airborne  during  sampling.  Additional  continuity  considerations 
Included  sampler  efficiencies  and  assay  techniques.  Barad  (1958a)  presented 


an  extensive  discussion  of  data  collection  and  assay  techniques  used.  He 
concluded  that,  for  the  S02  sampler  measurements,  there  was  a  cumulative 
error  of  about  ten  percent.  Sample  budget  calculations  from  the  ino-meter 
arc  vertical  array  also  indicate  very  efficient  tracer  recovery. 

Time-averaged  S02  concentration  measurements  taken  at  equally  spaced 
positions  along  each  arc  were  used  to  compute  CWID.  The  collected  concen¬ 
tration  data  were  averaged  over  a  ten-minute  period  for  the  data  report. 
Multiplying  the  reported  sampler  data  by  the  ten-minute  averaging  period 
produced  a  total  dosage.  The  sum  of  total  dosages  across  each  arc  multiplied 
by  sampler  spacing  Intervals  produced  CWIPs  which  were  then  used  in  equation 
(1)  to  obtain  values  of  sigma  z. 

Unfortunately,  cloud  centroid  measurements  were  not  made  at  each  arc 
during  Project  Prairie  Grass.  Except  for  the  100-meter  arc,  no  attempt  was 
made  to  measure  vertical  concentration  profiles.  Consequently,  the  behavior 
of  the  cloud  centroid  as  it  moved  beyond  100  meters  is  not  known.  To  solve 
equation  (1)  it  was  necessary  to  assume  no  cloud  centroid  rise  beyond  the 
100-meter  arc.  Although  necessary  to  obtain  solutions,  this  assumption  is 
not  well  founded.  Willis  and  Deardorff  (1976),  for  example,  have  documented 
rapidly  decreasing  surface  concentrations  in  unstable  conditions  due  to 
thermally  Induced  cloud  centroid  rise.  It  is  probable  that  cloud  centroid 
rise  affected  Project  Prairie  Grass  sampler  measurements,  especially  during 
unstable  conditions  at  the  more  distant  arcs. 

Meteorological  profile  data  accompanying  Project  Prairie  Grass 
diffusion  trials  were  collected  on  16-meter  towers,  each  instrumented  for 
wind  speed,  temperature,  and  wind  direction  measurements.  Also,  bi-vanes 
for  measurement  of  horizontal  and  vertical  wind  angles  were  mounted  on  a 
series  of  2-meter  masts.  These  data  were  collected  using  a  one-second 
averaging  time  for  twenty  minutes.  Data  processing  from  magnetic  tapes 
produced  temperature  gradients,  velocity  profiles,  and  wind  angle  standard 
deviation  information.  The  data  and  further  details  on  Project  Prairie  Grass 
are  in  Barad  (195Ra,b)  and  Haugen  (1959). 

V.  CHOOSING  A  TURBULENCE  PARAMETER 

Turbulence  data  accompanying  the  diffusion  trials  at  Project  Prairie 
Grass  consist  of  bi-vane  horizontal  and  vertical  wind  angle  measurements, 
presented  as  wind  angle  sigmas.  However,  the  use  of  vertical  wind  angle 
sigmas  derived  from  these  bi-vanes  is  not  widely  accepted.  The  primary 
objection  to  these  data  is  that  the  vanes  were  mounted  on  2-meter  masts. 
Because  of  proximity  to  the  ground,  it  would  appear  that  some  of  the 
turbulent  energy  would  be  truncated  and  the  data  would  therefore  not 
represent  the  full  amount  of  vertical  turbulence  acting  on  a  diffusing  cloud. 
This  truncation  would  be  most  noticeable  on  the  larger  eddies  generated  in 
the  unstable  regime. 

Given  the  limitations  of  the  bi-vane  sigmas  and  the  fact  that  only  16 
of  the  34  unstable  regime  trials  had  tabulated  sigmas,  an  alternative  means 
of  estimating  unstable  regime  vertical  turbulence  was  chosen.  This  led  to 
an  investigation  of  meteorological  profile  data  collected  during  the  Kansas 


1960  Field  Program  (Izumi,  1971).  This  project  was  carried  out  in  the  summer 
of  I960  at  a  field  site  in  southwest  Kansas.  The  Kansas  data  include  drag 
plate  measurements  for  u*  and  research  grade  wind  and  temperature  profile 
measurements  on  a  32-meter  tower.  Measurements  of  turbulent  velocity 
fluctuations  and  air  temperature  fluctuations  were  made  at  three  levels 
concurrently  with  profile  measurements.  The  basic  averaging  period  for  all 
data  was  15  minutes.  Turbulence  data  derived  from  these  measurements  do  not 
suffer  from  the  deficiencies  attributed  to  the  2-meter  bi-vane  data. 

A  relationship  was  established  between  the  standard  deviations  of  the 
vertical  velocity  (o  )  and  Kansas  unstable  regime  meteorological  profile  data 
using  the  following  assumptions: 

1.  A  log-linear  wind  velocity  profile  {equation  4)  describes  the 
horizontal  wind  speed  as  a  function  of  friction  velocity,  the  height  to 
roughness  logarithm,  and  a  diabatic  Influence  function  (ip).  The  unstable 


u  =  u*(ln(z/z0)  -  ^)/k 


(*) 


regime  diabatic  Influence  function  as  defined  by  Paulson  (1970)  is  presented 
in  equation  (5). 


ip  =  21n((l+*‘M/2)  +  ln((l+<f>"2)/2)  -  2tan"V  +  (5) 

2.  Empirical  equations  (6)  and  (7)  from  Dyer  (197^)  estimate  the 
unstable  regime  dimensionless  wind  shears  and  temperature  gradients. 

♦m  ■  (1  -  16(z/L))-*  (6) 

♦h  «  0  -  ISU/L))-5*  (7) 

Several  similar  equations  for  these  variables  are  found  in  the  literature. 
Equations  (6)  and  (7)  were  chosen  in  the  given  form  because  they  not  only 
fit  the  micro-meteorological  data  collected  at  Hay  and  Gurley,  New  South 
Hales,  Australia  (Dyer  and  Hicks,  1970),  and  provide  a  reasonable  fit  to  the 
unstable  regime  Kansas  data  as  corrected  by  Wieringa  (1980). 

Equations  (6)  and  (7)  imply  equality  between  R1  and  z/L  in  the  unstable 
regime.  Equations  (8)  and  (9)  thus  describe  the  relationships  between  the 
stability  parameters. 


R1  =  (z/L)4h/*2)  =  z/L  (8) 
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(9) 


a.  'a 


EU  =  1 00  ( z/ L )  /  ( (1  -  16(z/L))1'5(ln(z/z0)  -  *)) 


Equations  (6)  through  (9)  are  valid  for  a  z/L  range  of  0>z/L>-1.0.  Beyond  a 
z/L  of  -1.0,  Ri  becomes  steadily  larger  than  z/L.  A  sufficient  amount  of 
extremely  unstable  data  is  not  available  for  further  analysis  of  this 
phenomenon . 

With  the  relationships  given  by  equations  (4)  through  (6),  it  was 
possible  to  establish  a  link  between  Kansas  unstable  regime  profile  data 
and  sigma  w.  This  was  accomplished  by  examining  the  ratio  of  aw  to  u*.  This 
ratio  would  approach  unity  except  for  the  contribution  of  buoyancy  (gw?/e)  to 
vertical  turbulent  motions.  Using  Wieringa's  (1980)  suggested  adjustments 
to  the  Kansas  z/L  data,  it  was  found  that  the  ratio  of  o  to  u*  is  approx¬ 
imated  by  the  inverse  of  dimensionless  wind  shear.  Equation  (10)  also 
provides  a  close  fit  to  the  extremely  unstable  data  attributed  to  f-’onji 
(Businger,  1979). 


(10) 


For  an  absolutely  neutral  atmosphere,  equation  (10)  implies  equality  of 
o  and  u*.  However,  the  near  neutral  (0  >  z/L  >  -0.08)  measurements  at 
Kansas  indicate  a  to  u*  ratio  between  1.2  and  1.9,  which  is  consistent 
with  other  near  neutral  observations  (Pasquill,  1974).  The  situation  is 
further  clouded  by  Wieringa's  demurral  to  a  sigma  w  correction  factor  for 
the  Kansas  data.  A  a  to  u*  ratio  of  1.2  was  adopted  for  the  near  neutral 
condition,  pending  further  clarification. 

Using  equations  (4)  through  (10),  a  generalized  expression  is  found  for 
the  unstable  regime  by  linking  vertical  turbulence  (a  ,  radians)  to  meteor¬ 
ological  profiles.  Equation  (11)  describes  this  relationship.  For  near 


ae  =  aw/u  =  M-’VOnU/Zo)  -  *)  (11) 

neutral  conditions,  equation  (12)  pertains. 


oe  =  ow/u  =  1.2k/(ln(z/z0)  -  *)  (12) 


Mo  relationships  comparable  to  equations  (11)  and  (12)  were  found  for 
stable  regime  data.  However,  Waldron  (1980)  used  Project  Prairie  Brass  stable 
regime  data  to  establish  an  empirical  regression  equation  for  the  2-meter 
bulk  Richardson  number  and  sigma  e,  as  in  equation  (13). 


oe  =  EXP(-0 .2977BU  +  1  .416)  (13) 
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VI.  RESOLVING  VARIATIONS  IN  1 

Ideal  turbulence  measurements  would  be  made  by  a  perfectly  responding 
Instrument  during  a  steady  state  period  of  unbounded  random  turbulence. 

Cloud  growth  under  these  Ideal  conditions  would  be  described  as  a  linear 
function  of  the  travel  distance-turbulence  product.  In  reality.  Instrument 
response  Is  Imperfect  and  turbulence  varies  In  space  and  time.  To  compensate 
for  departures  from  Ideal  conditions,  equation  (14)  Introduces  a  variable 
correction  factor  (1).  Sigma  z  can  be  readily  computed  using  equation  (14) 


cz  =  aexi  (14) 

once  1  is  described  as  a  function  of  travel  distance  and  meteorological 
profile  data. 

To  depict  the  functional  relationship  of  1  to  travel  distance  and 
meteorological  profile  data,  a  field  of  1  values  distributed  over  a  range 
of  stabilities  and  travel  distances  was  needed.  This  field  of  1  values  was 
obtained  by  solving  equation  (14)  for  1  using  Project  Prairie  Grass  data. 
Sigma  z  and  sigma  e  values  were  computed  at  each  trial  arc  distance  using 
procedures  described  In  III  and  V  above.  The  correction  factor  computed  In 
this  way  compensates  for  departures  from  Ideal  turbulence  measurements  and 
deviations  from  the  assumptions  implied  by  equation  (1).  This  procedure  is 
useful  so  long  as  the  resultant  field  of  1  "alues  can  be  described  as  a 
function  of  the  parameters  responsible  for  producing  the  measured  dosages. 

It  was  assumed  that,  with  corrections  for  continuity,  source  configuration, 
and  roughness,  the  relationships  derived  using  Project  Prairie  Grass  data 
would  be  applicable  to  other  data  sets. 

Real  values  of  1  were  calculated  for  all  arc  distances  on  Project 
Prairie  Grass  trials  for  which  data  were  available.  A  derivation  In  Appendix 
A  shows  that  the  upper  limit  of  real  i  values  is  reached  when  CWID*u/Q  equals 
0.323.  Beyond  this  limit  values  of  1  are  imaginary  and  cannot  be  used  for 
evaluation. 

Analysis  of  the  field  of  1  values  Indicated  a  systematic  variation  with 
travel  distance,  but  with  considerable  scatter  between  trials,  ^uch  of  this 
scatter  was  found  to  be  a  function  of  a  non-dimensional  stability  parameter, 
the  bulk  Richardson  number.  The  2  meter  values  of  BU  provided  the  best  fit. 
Two  meter  bulk  Richardson  numbers  were  computed  from  profile  data  using 
equation  (15). 


BU  =  100gz2(AT/AZ  +  0.01 )/Tu2 


(15) 
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where: 


T  =  2  meter  temperature  (degrees  Kelvin) 

u  =  2  meter  wind  speed  (meters  per  second) 

z  =  height  of  wind  speed  measurement,  2  meters 

aT/az  =  %  to  4  meter  temperature  differential  (degrees  Celsius  per  meter) 

Puring  analysis  of  the  relationship  between  i  and  BU  it  became  apparent 
that  the  variation  of  i  with  RU  assumed  an  arctangent  shape.  The  arctangent 
shape  allows  for  large  changes  in  i  as  Btl  passes  through  near  neutral 
stability,  with  decreasing  changes  in  i  as  PIJ  approaches  the  stability 
extremes.  Also,  the  amplitude  of  the  fitted  arctangent  curves  appear  to 
increase  systematically  with  distance.  Equation  (16)  was  developed  to 
describe  this  expanding  topological  field  of  1  values  over  a  continuum  of 
stability  and  distance. 


1  =  (-ATAN(a*BU  +  b)  +  c)/d  (16) 

where: 

a  «  10.0  (17) 

b  =  x0*21  -  0.9  (lo) 

c  «  3900x'1*6  +  1  .5705  (19) 

d  =  1860X'1-2  +  0.37  (20) 


Variables  a  through  d  in  equation  (16)  allow  for  adjustment  o*  the 
expanding  arctangent  field.  Equations  (17)  through  (20)  define  a  "best  fit" 
field,  subject  to  several  limiting  conditions.  The  first  limiting  condition 
is  that  the  coefficient  "a"  must  be  constant  for  all  distances.  Any  variation 
in  "a”  with  distance  would  imply  a  cross  correlation,  or  a  functional  relation¬ 
ship,  between  Oil  and  travel  distance.  The  second  limiting  condition  is  that 
i  can  never  become  zero  or  negative.  Because  an  arctangent  is  limited  in 
range  of  radian  values  to  ±1.5708,  the  minimum  value  of  "c"  must  be  1.5708. 

The  remaining  coefficient  values  were  selected  by  trial  and  error  to  provide 
the  best  fit  to  the  i  values  derived  from  Project  Prairie  Grass  data. 

Figure  1  depicts  a  segment  of  the  i  field  given  by  equation  (16)  for  arc 
distances  of  50  to  800  meters  and  a  BU  range  of  ±2.0.  The  most  noticeable 
feature  on  this  figure  is  the  dramatic  change  in  i  for  small  changes  in  the 
near  neutral  values  of  BU,  illustrating  the  marked  transition  in  diffusing 
power  as  the  atmosphere  passes  between  stable  and  unstable  regimes.  Figure  1 
also  illustrates  expansion  of  the  arctangent  relationship  between  BU  and  i  as 
travel  distance  Increases.  Also,  Figure  1  shows  correction  factor  values  to 
be  small  and  decreasing  with  travel  distance  for  the  stable  regime,  and  large 
and  increasing  with  travel  distance  for  the  unstable  regime.  Equations  (17) 
through  (20)  allow  for  infinite  adjustment  of  this  topological  field  to  fit 
the  array  of  Project  Prairie  Grass  1  values. 
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VII.  VERIFICATION  ON  DEPENDENT  DATA 

As  a  check  on  the  sigma  z  computation  procedure,  a  verification  was 
performed  with  the  Project  Prairie  Grass  data.  Sigma  z  values  were  computed 
from  equation  (14)  using  profile  derived  sigma  e  values  and  correction  factors 
from  equation  (16).  These  sigmas  were  then  used  In  equation  (1)  for  re¬ 
computation  of  crosswind  Integrated  dosages,  designated  as  CWIP_.  The  CWID 
were  then  compared  to  the  original  CWIDs  obtained  from  sampler  arc  average  p 
concentration  measurements.  Fractional  error  (FE)  defined  by  equation  (21) 
was  used  as  the  basis  for  comparison  of  CWIDp  to  CWID. 

FE  =  (CWIDp  -  CKID)/Jj(CWIDp  +  CHID)  (21) 

Fractional  error  Is  the  difference  between  observed  and  predicted  CWIDs, 
normalized  by  the  average  of  observed  plus  predicted  CWIDs.  Horst  (1979) 
described  FE  computed  in  this  way  as  providing  a  logarithmically  unbiased 
means  of  presenting  the  fractional  difference  between  observed  and  predicted 
parameter  values.  The  maximum  possible  fractional  difference  between  two 
parameters  produces  a  FE  of  2.D,  and  a  factor  of  2  difference  is  indicated 
by  a  FE  of  0.667.  Mean  fractional  error  (FE)  is  defined  by  equation  (22), 
where  N  is  the  number  of  trials  evaluated.  Root  mean  square  fractional  error 


FE  =  £FE/N 


(22) 


(FErms)  is  defined  by  equation  (23). 

FErms  =  (JFE^N)*5  (23) 


There  were  71  Project  Prairie  Grass  trials.  No  S02  concentration  data 
were  available  for  trials  48S2,  63,  and  64.  With  five  arc  distances  per  trial, 
340  cases  remained.  Concentration  data  were  missing  at  800  meters  and  400 
meters  for  trials  39  and  40  respectively,  and  trial  3  was  eliminated  because 
most  of  the  cloud  missed  the  samplers,  leaving  a  total  of  333  cases  for 
evaluation.  Results  of  the  fractional  error  evaluation  for  the  remaining 
cases  are  presented  in  Table  1.  Pasquill  stability  categories  used  in  Table  1 
were  determined  by  Colder 's  (1972)  procedure.  Results  presented  in  Table  1 
compare  favorably  with  those  obtained  by  Horst  (1979)  in  his  dimensionless 
parameter  treatment  of  the  same  data.  Less  than  4  percent  of  the  CWIDp 
differed  from  CWID  by  a  factor  of  2  or  more.  The  largest  fractional  errors 
were  generally  found  at  the  800-meter  arc  and  stability  categories  A  and  F. 
Overpredictions  of  CWIDp  during  unstable  conditions  at  800  meters  were 
probably  due  to  cloud  centroid  rise.  Overpredictions  of  CWIPp  durinq  stable 
conditions  may  be  due  to  cloud  travel  along  a  laminar  sublayer,  skirting  the 
samplers  at  1 .5  meters . 
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Table  1.  Mean  and  Root  Mean  Square  Fractional  Error  for  Project  Prairie 
Grass  Data  at  Indicated  Distances  and  Pasquil1  Stability  Categories. 


Distance 
(meters ) 

Mean 

Fractional  Error 

Pasquill  Stability  Category 

A 

B 

C 

B 

E 

F 

ALL 

50 

-0.17 

-0.08 

-0.15 

-0.01 

-0.12 

0.28 

-0.03 

100 

0.02 

-0.06 

-0.07 

0.15 

-0.04 

0.15 

0.05 

200 

0.07 

-0.26 

-0.06 

0.18 

0.03 

0.11 

0.06 

400 

-0.04 

-0.25 

0.01 

0.11 

-0.09 

0.07 

0.02 

800 

0.40 

-0.05 

0.15 

-0.13 

-0.40 

0.15 

-0.03 

ALL 

0.05 

-0.14 

-0.02 

0.06 

-0.12 

0.15 

0.01 

Root  Mean 

Square 

Fractional 

Error 

Distance 

(meters) 

Pasquill  Stability  Category 

A 

B 

C 

D 

F 

ALL 

50 

0.37 

0.47 

0.22 

0.09 

0.13 

0.42 

0.26 

100 

0.36 

0.43 

0.20 

0.19 

0.12 

0.33 

0.25 

200 

0.35 

0.43 

0.27 

0.23 

0.12 

0.36 

0.28 

400 

0.26 

0.47 

0.32 

0.19 

0.18 

0.55 

0.32 

800 

0.57 

0.65 

0.29 

0.27 

0.46 

0.45 

0.39 

ALL 

0.39 

0.50 

0.26 

0.20 

0.24 

0.42 

0.30 

1? 


VIII.  VERIFICATION  ON  INPEPENPENT  DATA 

After  successful  application  to  Project  Prairie  Grass  data,  the  analytic 
method  was  applied  to  independent  sets  of  data.  The  first  of  these  consisted 
of  twenty  three  surface  source  trials  conducted  at  the  Horizontal  Grid  on  the 
U.  S.  Army  Pugway  Proving  Ground,  Utah.  These  trials  are  known  as  the  PC  41? 
Trials  (Cramer  et  aJL_,  1965).  A  major  advantage  of  using  BC  412  data  for 
verification  is  that  a  large  number  of  vertical  concentration  profiles  were 
obtained  from  an  extensive  vertical  array  of  samplers.  Pirect  vertical 
concentration  measurement  obviated  source  strength  considerations,  as  the 
sigma  z  was  determined  directly  from  the  vertical  profiles,  reposition  was 
assumed  to  be  insignificant  for  the  short  travel  distances,  as  1  -  10  ym 
Bacillus  subtilis  var  niger  particles  were  used  as  a  tracer.  It  was  also 
assumed  that  al 1  samplers  had  similar  efficiencies  and  were  subjected  to 
similar  handling  and  assay  techniques.  Source  generation  was  accomplished 
using  trailer  mounted  disseminators  pulled  along  a  track  perpendicular  to  the 
wind  direction.  It  was  determined  that  this  dissemination  procedure  produced 
a  cloud  with  a  source  sigma  of  1.75  meters  and  a  release  height  of  approx¬ 
imately  1.5  meters  above  ground  level. 

Meteorological  profile  data  accompanying  the  PC  412  Trials  were  collected 
on  30  meter  towers  instrumented  for  temperature  and  wind  profile  measurements. 
The  h  to  4  meter  temperature  differentials  were  determined  from  trial  profile 
data,  and  2  meter  bulk  Richardson  numbers  were  calculated.  The  unstable 
regime  BU  were  processed  through  equations  discussed  in  Section  V  to  produce 
estimates  of  sigma  e.  Estimates  of  sigma  e  for  the  stable  regime  trials 
were  obtained  using  equation  (13)  modified  by  a  roughness  correction  (Smith, 
1972).  The  roughness  correction  is  the  ratio  of  PPG  Horizontal  Grid  to 
Project  Prairie  Grass  z/z0  logarithms.  For  a  PPG  Horizontal  Grid  z0  of  2.3 
centimeters,  a  Project  Prairie  Grass  z0  of  0.9  centimeters,  and  a  height  of 
2  meters,  the  roughness  correction  is  given  by  equation  (24).  The  stable 


ln(2/0.023)/ln (2/0.008)  =  1 .24 


(2«) 


regime  sigma  e  regression  equation  with  the  roughness  correction  are  given 
by  equation  (25).  Sigma  z  was  computed  as  the  square  root  of  the  sum,  of 


ae  =  (FXP(-C.2977BU  +  1 .41 6) ) -1 .24  (25) 

the  source  and  turbulent  growth  variances,  as  provided  In  equation  (26). 

az  =  (1.752  +  (OgXi)2)*5  (26) 
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Results  of  verification  with  BC  412  data  are  presented  in  Table  2. 
Fractional  errors  were  computed  for  the  predicted  versus  measured  sigma  z 
using  the  procedure  applied  during  Project  Prairie  Grass  CWIPp  verification. 
The  FErms  for  the  entire  data  set  was  0.31,  and  6  percent  ofpthe  predicted 
sigma  z  differed  from  measured  values  by  more  than  a  factor  of  2. 

Table  2.  Mean  and  Root  Mean  Square  Fractional  Errors  for  Computed  versus 
Observed  Sigma  Z,  DPG  Horizontal  Grid  BC  412  Data. 


Distance 

(meters) 

Unstable  Regime 

FE  FErms 

N 

Stable  Peqime 

FE  FErms" 

N 

30 

-0.20 

0.25 

g 

-0.09 

0.32 

15 

91 

-0.05 

0.15' 

8 

-0.19 

0.32 

15 

274 

0.36 

0.55 

5 

-0.13 

0.28 

15 

427 

0.22 

0.34 

3 

0.05 

0.26 

15 

793 

— 

— 

0 

0.25 

0.35 

14 

ALL 

0.02 

0.33 

24 

-0.O3 

0.31 

74 

A  second  evaluation  on  an  independent  data  set  was  performed  using  data 
collected  at  the  National  Reactor  Testing  Station  (NRTS)  in  southeast  Idaho 
(Islitzer  and  Dumbauld,  1963).  Diffusion  measurements  were  made  with  sampling 
arcs  operated  100,  200,  400,  R00,  1600,  and  3200  meters  downwind  from  the 
source.  The  source  was  a  uranine  dye  in  small-particulate  form  released  at  a 
height  of  1  meter  for  60  minutes  per  trial.  Crosswind  integrated  concentration 
data  were  reported  for  each  arc,  normalized  to  a  source  strength  of  1  gram 
per  second  and  a  transport  wind  of  5  meters  per  second.  The  cloud  centroid 
and  measurement  height  were  both  assumed  to  be  at  the  surface  (z  =  h  =  0.0). 

The  NRTS  test  site  was  described  as  flat  with  widely  distributed  sage¬ 
brush.  Roughness  was  estimated  at  1.5  centimeters.  Meteorological  profile 
measurements  of  temperature  gradients  and  wind  velocities  were  made  near  the 
tracer  release  point.  Tabulated  4-and  8-meter  Ri  were  converted  to  z/L  and 
reduced  to  2-meter  values.  For  the  unstable  regime,  this  allowed  the  use  of 
the  same  stability  parameters  as  were  applied  to  Project  Prairie  Grass  data. 

No  stable  regime  evaluations  were  performed  because  insufficient  data  were 
available  for  calculation  of  comparable  stability  parameters. 

Evaluation  of  continuity  considerations  for  the  NRTS  trials  was  difficult 
because  tracer  deposition  was  suspected  (Islitzer  and  Dumbauld,  1963). 

Equation  (27),  attributed  by  Islitzer  and  Dumbauld  to  Chamberlain,  was  used 
to  compensate  for  loss  of  source  strength  with  downwind  travel.  Islitzer  and 
Dumbauld  found  a  V0  of  5.2  centimeters  per  second  and  a  nz  of  -0.24  to  be 
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qc  =  Q.EXP(-4VQx/((/27r)nzuaz)) 


(27) 


where: 

Qc  =  adjusted  source  strength  (grams) 

Vo  =  deposition  velocity  (meters  per  second) 

n2  =  dispersion  parameter  (dimensionless) 

average  values  for  these  parameters.  Sufficient  data  were  not  available  to 
use  equation  (3)  for  u  calculation.  Therefore,  the  Islitzer  and  numbauld 
tabulated  wind  speed  data  were  used  for  u  In  deposition  calculations. 

Table  3  provides  the  results  of  fractional  error  evaluation  for  observed 
versus  predicted  crosswind  Integrated  concentrations  with  unstable  regime 
NRTS  data. 

Table  3.  Mean  and  Root  Mean  Square  Fractional  Errors  for  Computed  versus 
Observed  Crosswind  Integrated  Concentrations,  NRTS  Pata. 


Distance 

Unstable  Regime 

(meters ) 

FE 

FErms 

N 

100 

-0.03 

0.45 

14 

200 

-0.26 

0.48 

23 

400 

-0.16 

0.48 

13 

800 

-0.05 

0.45 

22 

1600 

0.13 

0.45 

14 

3200 

n.4P. 

0.68 

14 

ALL 

-0.C1 

0.50 

100 

The  FErms  of  0.50  Is  higher  than  that  for  the  Project  Prairie  Grass  or 
the  BC  412  data,  but  Is  well  within  a  factor  of  two.  Twenty  percent  of  the 
cases  were  In  error  by  a  factor  of  two  or  more.  The  largest  errors  were 
found  at  the  3200-meter  distance,  where  systematic  overprediction  occurred. 
These  overpredictions  may  be  due  to  insufficient  deposition  calculated  with 
equation  (26),  insufficient  cloud  growth  prediction  using  i,  or  some  combin¬ 
ation  of  the  two. 
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IX.  CONCLUSIONS 

An  analytic  method  for  calculating  sigma  z  from  meteorological  profile 
data  has  been  developed  and  demonstrated.  With  appropriate  adjustments  for 
source  sigmas,  deposition,  and  surface  roughness,  the  method  was  successfully 
used  on  Independent  data  sets,  consistently  providing  estimates  of  sigma  z 
or  crosswind  integrated  dosage  well  within  a  factor  of  2  accuracy. 

The  sigma  z  calculation  method  is  based  on  the  use  of  a  correction 
factor  adjustment  to  rectilinear  cloud  growth.  For  Project  Prairie  Grass 
data,  this  correction  factor  was  a  function  of  stability  and  travel  distance, 
described  over  a  continuum  of  stability  and  travel  distance  by  an  expandina 
arctangent  equation.  Because  cloud  growth  at  other  locations  observes  the 
same  functional  relationship  to  stability  and  travel  distance,  the  method 
has  universal  applications.  One  must  note,  however,  that  the  correction 
factor  field  developed  in  this  report  is  based  on  a  surface  source  and  sampling 
at  1.5  meters  above  ground  level.  Changes  in  source  or  sampling  height  may 
require  adjustments  to  the  correction  factor  field  pattern.  Sufficient  data 
to  define  such  adjustments  are  not  available. 

The  sigma  z  calculation  method  requires  an  estimate  of  vertical 
turbulence,  and  procedures  for  estimating  vertical  turbulence  from  meteor¬ 
ological  profile  data  are  demonstrated.  The  turbulence  estimation  procedure 
for  the  unstable  regime  is  valid  through  a  z/L  of  -1.0.  Application  of  this 
procedure  in  extremely  unstable  conditions  where  z/L  exceeds  -1.0  is 
uncertain  because  adequate  relationships  between  stability  parameters  have 
not  been  developed  for  these  conditions.  Likewise,  use  of  the  regression 
equation  which  describes  sigma  e  as  a  function  of  Bl.1  becomes  uncertain  beyond 
a  BU  of  +4.0.  Also,  relationships  between  stability  parameters  are  not  well 
established  for  the  stable  regime,  requiring  a  known  2-meter  Bl!  before  stable 
regime  sigma  e  values  can  be  estimated. 

Although  the  sigma  z  calculation  method  is  applicable  over  the  full 
range  of  stability  and  travel  distance,  present  procedures  for  estimating 
vertical  turbulence  are  subject  to  revision.  Fortunately,  coefficients  of 
the  arctangent  equation  are  infinitely  adjustable.  As  new  turbulence 
computation  procedures  are  developed,  coefficients  of  the  arctangent  equation 
will  allow  for  adjustments  to  accomodate  revisions  while  retaining  the 
functional  form. 
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APPENDIX  A. 


DERIVATION  OF  A  CWIDU  LIMIT  FOr  THE  CWID  EQUATION 


The  basic  Cwid  equation  including  reflection  terms  is 


C“ID  =  jAr  + 

/2ttozU  20^  z 


Substituting  owti  for  az  and  rearrangement  of  terms  yields: 


CWID  0wtu  _  if  /-(z-h)2  „f-(z+h)2n 

- q - 


(Al) 


(A2) 


with 


(z+h)2  =  (z-h)2  +  Uzh 

and  defining 


A  = 

B  = 


1 

2ow2t 2  * 

✓57  CWID  gvtU 


Q 


then. 


(A3) 


(A4) 

(A5) 


B  =  1  [e(-AL2-h.Ll  +  e(~AK-z-h)itiizh))]  (A6) 

i  i2  i2 

Considering  i  a  general  differentiable  variable,  and  simplifying  equation 
(A6)  by  defining 

C  =  A(z-h)2,  (A7) 

D  =  a( ( z-h )2  +  Uzh) ,  (A8) 

F(i)  =  i-  [e^_q/l2^  +  e^-D/i2^]  -B  =  0  (A9) 

1 


A-l 


Taking  the  derivative. 


df(i)  =  d_(l)[e(-c/i2)  +  e(-D/i2 )  -j  +  l_  d_^(-C/i2)  + 
di  di  i  i  di 


At  the  maximum  value  of  a  function,  the  derivative  is  zero 

df ( i )  =  o 

di 

Simplifying  equation  A10  and  gathering  terms, 

i2  df ( i )  B  (2C  -  i)  e(“C/i2)  +  -  !)  e(-D/i2 )  =  0 

di  i2  i 


For  a  sampling  height  (2)  of  1.5  meters  and  cloud  centroid 
0.5  meters, 

(z-h)2  =  1.0, 

Uzh  =  3.0 

Then, 


df(i) 

di 


(-UA/i2) 


=  0 


As  long  as  NA/i2  is  small, 
em/i*  .  1  *  NA/i2 


Then, 

(2A/i2  -  l)(l  -  A/i2)  +  (8A/i2  -  1 ) (l  -  UA/i2)  =  0 


Simplifying, 

-3hk2  +  15Ai2  -  2i4  =  0 
Solving  equation  (A18)  quadratically , 
i2  .  1|A  ±  SJT 


A- 2 


e(“D/i2)] 

(A10) 

(All) 

(A12) 

height  (h)  at 

(A13) 

(A14) 

(A15) 

(A16) 

(A17) 

(A18) 

(A19) 


It  can  be  shown  by  De  Moivre's  theorum  that 


i4  =  (1U.0625  +  2.93T)A2  (A20) 

Hence: 

i2  =  1».123A  (A21) 

Substituting  equation  (A21)  into  equation  (A6),  along  with  equations 
(A13)  and  (A14)  yields 

B  =  0-5731 

StT  (A22) 

From  equation  (A4)  one  obtains 

-i  =  ^owt  (A23) 

Then,  employing  equations  (A22),  (A23),  and  (A5), 

B  =  Q.-JI.31  =  o.5731^cr  t  =  qwtu 

Jk  Q  (A24) 

Or, 

cyiD  y  =  0.3233,  (A25) 

Q 


which  is  very  close  to  the  exact  solution  achieved  iteratively  (0.3235)- 

Any  CWID  u/Q  value  >0.3233  violates  the  limit  [equation  (All)]  on 
equation  (A10)  and  therefore  provides  no  real  solution  to  the  CWID 
equation.  The  limit  presented  in  equation  (A25)  is  unique  for  the 
given  values  of  z  and  h,  but  the  procedure  can  be  used  to  solve  for 
other  z  and  h  values.  For  a  condition  where  the  derivative  is  nonzero, 
it  is  apparent  that  there  are  two  possible  solutions  for  equation  (A19). 
This  situation  can  be  resolved  by  choosing  the  smaller  real  value  for 


